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Abstract 

We present a class of realistic unified models based on supersymmetric 5*0 (10) wherein 
issues related to natural doublet-triplet (DT) splitting are fully resolved. Using a minimal set 
of low dimensional Higgs helds which includes a single adjoint, we show that the Dimopoulos- 
Wilzcek mechanism for DT splitting can be made stable in the presence of all higher order 
operators without having pseudo-Goldstone bosons and flat directions. The /i term of order 
TeV is found to be naturally induced. A Z2-assisted anomalous U{1)a gauge symmetry plays 
a crucial role in achieving these results. The threshold corrections to a3{Mz), somewhat 
surprisingly, are found to be controlled by only a few effective parameters. This leads to a 
very predictive scenario for proton decay. As a novel feature, we find an interesting correlation 
between the d = 6 {p ^ e^ir^) and d = 5 (p — t- VK^) decay amplitudes which allows us to 
derive a constrained upper limit on the inverse rate of the e'^'K^ mode. Our results show that 
both modes should be observed with an improvement in the current sensitivity by about a 
factor of five to ten. 
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1 Introduction 



Although yet to be seen, proton decay is an indispensable tool to probe nature at truly high energies 
(~ 10^^ GeV). It still remains as the missing piece of grand unification [1, 2, 3]. In the hght of 
a new set of planned detectors including those at the forthcoming deep underground laboratory 
DUSEL [4] and the HyperKamiokande, we propose to address in this paper certain well known but 
partially unresolved theoretical issues of supersymmetric (SUSY) grand unification (GUT) which 
are especially relevant to proton decay. 

Strong empirical support for grand unification arises not only from the observed quantum num- 
bers of quarks and leptons and the quantization of electric charge, but in particular from the meeting 
of the three gauge couphngs at a scale Mgut ~ 2-10^^ GeV that occurs in the context of low energy 
SUSY [5], and the tiny neutrino masses, as observed in neutrino oscillation experiments. The latter 
fit extremely well with GUT symmetries that include the symmetry SU{2)l x SU{2)r x SU{4)c 
[1], the minimal such symmetry being SO{10) [6]. We will therefore discuss proton decay in the 
context of supersymmetric SO{10). The purpose of the present paper is to pay special attention 
to the problem of the so-called doublet-triplet (DT) splitting and to study the imphcations of its 
resolution for proton decay. 

The DT splitting problem is common to all grand unified theories based on simple gauge groups. 
In SUSY SO{10) models the two Higgs doublets of MSSM, a color triplet and an anti-triplet lie 
(typically) in a 10-dimensional representation H{10). The color triplets need to be superheavy so as 
to avoid rapid proton decay and also to preserve gauge coupling unification. Keeping the doublets 
light and the triplets superheavy self-consistently is the doublet-triplet splitting problem. 

A natural solution to this problem, avoiding severe fine-tuning is realized in SUSY 5*0(10) by 
the so called Dimopoulos-Wilczek (or the missing VEV) mechanism [7]. It involves a coupling of 
two 10-plets of the form H (10) A{45)H' (10) with the adjoint A(45) having a GUT scale VEV in 
the {B — L)-preserving direction: 

{A) = i(72 <H) Diag (a, a, a, 0, 0) . (1) 

This structure contributes to the triplet and not to the doublet masses, and thereby can lead to 
natural DT splitting without fine-tuning. 

Given the very large hierarchy between the doublet and triplet masses, however, one must 
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ensure: (i) that the missing VEV pattern for A{45) in Eq. (1) is stable to a high enough accuracy 
in the presence of all allowed higher dimensional operators; (ii) that there are no undesirable pseudo- 
Goldstone bosons; and (iii) that there are no flat directions which would lead to VEVs of flelds 
undetermined. Furthermore, (iv) one must also examine, by including all GUT-scale threshold 
corrections to the gauge couplings, the implication of the doublet-triplet splitting on coupling 
unification and on proton decay. To our knowledge, while some of these issues have been partially 
addressed in the literature (e.g. see [8], [9], [10]), and major progress was made in Ref. [11] with 
regard to the issues (i) and (ii), simultaneous resolution of all four issues has so far remained a 
challenge. 

In this paper we present a predictive class of 5*0(10) models, based on a minimal Higgs system, 

in which all the issues of DT splitting mentioned above are resolved, and where the threshold 

corrections to the gauge couplings and their implications for proton decay are properly studied 

as well. The Higgs sector we consider has a single adjoint, along with vectors and spinors. Such 

a low dimensional Higgs system would lead to smaller threshold effects [8,10], unlike in models 

[8, 10] which employ multiple adjoints and/or 54 dimensional Higgs. A postulated Z2-assisted 

anomalous U{1)a symmetry (which may have a string origin [13]) plays a crucial role in obtaining 

our results. We find somewhat surprisingly that the GUT scale threshold corrections to a3{Mz) are 

determined in terms of a very few parameters. This makes the model rather predictive for proton 

decay. As a novel feature, we find an intriguing correlation between the d = 6 and d = 5 proton 

decays, which respectively lead to p ^ e+vr" and p VK^ as the dominant decay modes. The 

correlation is such that the empirical lower limit on r~^(p VK^) provides a constrained upper 

limit on V~^{p — > e+7r°). Our results show that both decay modes should in fact be discovered with 

an improvement in the current limits on lifetimes by about a factor of five to ten. 

^An alternative class of 50(10) models utilizing larger dimensional (e.g. 126) Higgs fields has been studied in 
Ref. [12]. These models have the interesting feature that -R-parity is automatic, being part of the gauge symmetry. 
However, threshold corrections are rather large in these models, making quantitative predictions for az{mz) and 
proton decay difficult (see attempts in this regard by Aulukh and Garg [12]). 
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2 Stabilizing doublet-triplet splitting 



In order to break SUSY SO{10) to the supersymmetric standard model with a stabihzed DT 
sector, and for the subsequent breaking of the electro-weak symmetry, we shall use a minimal 
low dimensional Higgs system. It consists of a single adjoint A(45), two pairs of spinor-antispinor 
superfields {C(16) + (7(16)} and {C"(16) + C"(T6)}, two 10-plets H{10) and H'{10), as well as two 
SO{10) singlets S and Z. The second spinorial pair C + C' is introduced, following Ref. [11], 
to avoid pseudo-Goldstonc degrees of freedom while maintaining the Dimopoulos-Wilczek VEV 
structure for A (cf: Eq. (1)). The S and Z superfields are needed to fix various VEVs in the 
required directions through their superpotential couplings. 

We supplement the gauge symmetry by a Z2-assisted anomalous U{1)a symmetry in order to 
stabihze the VEV pattern of Eq. (1) [14], [10]. The charges of the Higgs fields and those of the 
three matter families 16j under U{X)a x Z2 are listed in Table 1. Here /c is a positive integer which 



Table 1: IA{\)a and Z2 charges Qi and Ui of the superfield 0j. The transformations under IA{1)a 
and Z2 are respectively 0j e*'^'0j and 0j e'^'^^^j. 





A(45) 


H{10) 


H'{10) 


C{16) 


C{16) 


z 


s 


C"(16) 


C'{T6) 


16l,2 


I63 


Q 





1 


-1 


2k 


1 
2 


2 
k 


2 
k 


2k 


k+8 
2k 


Ql,2 


1 
2 


UJ 


1 





1 








1 











A,2 






is unspecified for the moment. The superpotential of the symmetry breaking sector, consistent with 
these symmetries, isW — Wi + W2 + W3, where 

W, = MAtrA^ + ^ (trAf + ^trA^ , (2) 

W2 = c(^ZA + ^CC + c,s]c' + C'i^ZA + ^CC + C2s\c, (3) 

= \iHAH'+(\h'SZ''-^ + X'h^Z'')^^ + \2HCC+^AH'CC' . (4) 

For simphcity we assume that the S'O(IO) contractions in the CC terms with coefficients 61,2 in Eq. 
(3) are in the singlet channel. In the second term of Eq. (4) the operator Z^ can appear only when 
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k is even. While our consideration of DT splitting will hold for all k, if k is odd, matter parity 
is automatic, being part of U{1)a- The choice of k — 5, which we will use, is phenomenologically 
preferred, in particular for suppressing adequately all d — 5 proton decay operators including those 
induced by Planck scale physics. Higher order operators such as A^/M^ etc. are not exhibited in 
Eqs. (2)- (4) because they are inconsequential for our purposes. The charges gi.2 and the parities 
Pi^2 of the first two families are left unspecified for the present. They will however be relevant for 
the generation of quark and lepton masses. 

Typically, we expect that the non-renormalizable operators such as and A^-terms would be 
induced by quantum gravity effects involving exchange of heavy states in the string tower. Thus, we 
expect the cut-off scale ~ Mp\ or Mstj-mg ~ 10^*^ GeV. We shall take all dimensionless couplings 
to be of order unity, i.e., in the range (1/4 — 2). 



Using the SUSY preserving condition Fz — Fg — 0, together with the choice (C) — (C 



c, (A) 7^ (which is one allowed option among the discrete set of degenerate vacuum solutions) , we 
get (CC") = {CC) = and (C) = (C") = 0. The VEV of A is then determined entirely by Wi of 
Eq. (2). Setting Fa = 0, we find a solution in the B — L direction as in Eq. (1), with 



With Aa, A'^ ~ 1 and ~ 10^^ GeV, we need to choose Ma ~ 10^^ GeV to obtain a ~ Mqut ~ 
2 ■ 10^® GeV. Demanding F-flatness conditions Fq' = Fq, = and using the notations z = (Z) and 
s = (S) we get s = 4-pi and z = ^p2, where pi = ^^"^-^^"^ , 02 = ^'"■'-^^"^ . We note that for all 

\ ' ° M^r'^ 3a"^' aiC2— a2Ci ' aiC2—a2Ci 

dimensionless couplings in the Lagrangian being in the range (1/4 — 2), the effective couplings pi,2 
can naturally take values as small as about 1/50. 

The sum of the VEVs gets further constrained as follows. The anomalous U{1)a symmetry, 
presumed to have a string origin, generates the Fayct-lliopoulos term ^ through quantum gravity, 
which is given by [13] ^ = ^192^2' TrQ^, where Qst denotes the string coupling and Mpi ~ 2.4-10^* GeV 
is the reduced Planck mass. In our model, the particle spectrum of Table 1 would lead to Tr((5^) = 
-8 - 60/A; + 16(gi + ^2) = -84/5, (for A; = 5, qi,2 = -1/2 + 3/k, see later). This value will 
however be modified if there are additional singlets in the full theory. (Semi-realistic string solutions 
[15], possessing an anomalous U{1)a, typically lead to |Tr(5A| ~ 30 — 100.) With the charges 
in Table 1, the vanishing of Da = ^ + J2iQi\{4>i)\'^ = (required for preserving SUSY), yields 




(5) 



2{6Xa + X'a) ■ 
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+ |zp + |sp = — I'C- Thus, the VEVs of all the fields get determined. We see that quite naturally, 
the VEVs c, 2; ~ (few — 10) x Mqut, and s ~ (10~^ — 10~^) x Mgut can arise, with the precise values 
depending on the order one couplings. Let us note that this setup also allows for additional singlet 
fields {Pi} which can play a role in the Da = condition (for Pj with positive U{1)a charges) and 
can modify these estimates somewhat, without upsetting the stability of DT splitting. 

Substituting the VEVs of the heavy fields in Eqs. (3) and (4), we derive the mass matrices 
and Mt for the SU{2)l doublets and 5'C/(3)c-color triplets (written in the SU{5) notation): 



5c 



5c' V 



5h 


5h' 


5c 


5c' 





r]D,T>^ia 


A2C 





-r]D,T>^lO' 


Mh' 

















K'D,tYi 





Yd,t 







(6) 



with {vd,Vt) = (0,1), {kd,kt) = (3,2). Here Mh' ^ (Xh'SZ''-^ + X'H'Z')/Mt' (see Eq. (4)), 
=2ai.2-2a/(M*) and Yd,t ~ A3(^)c/(M*). For k = 5 the dominant contribution to Mh' comes 
from the operator Xh'Sz'^/M^, which is in the range (10^^ — 10^^) GeV. The suppressed mass of 
H' will be crucial for an adequate suppression oi d — 5 proton decay. The entry Mc in Eq. (6) 
(allowed by the stability of Higgs doublet mass) would arise from the operator SZ'^C'C' /M"^ and 
yields Mc ~ (10~^ to 10~^) x Mqut if M ~ which happens if the superfields that are integrated 
out have GUT scale masses. 

The zeros in the first column of Eq. (6) are ensured, in the presence of all higher dimensional 
operators, for the doublet mass matrix by the U{1)a x Z2 symmetry. The main reason for this 
all-order stability of the Higgs doublet masses is that all the effective Higgs fields (i.e. any positive 
power of Z, 5" and CC) which have super-large VEVs are positively charged under W(1)a, and 
can not couple to if^, which is also positively charged. Thus, with r^D = one pair of the Higgs 
doublets will be massless, while the remaining three pairs of doublets become superheavy. The role 
of the Z2 symmetry is that it allows the coupling of H to H' only through A (or odd powers of A). 
Such couphngs, however, do not generate a doublet mass due to the VEV structure in Eq. (1) of 
{A). The VEV pattern of {A) along the B — L direction is also guaranteed to be stable because 
of the U{1)a symmetry. Indeed, note that the symmetry U{1)a does not allow any superpotential 
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coupling involving A, C and C of the form A"(CC')'". It is only these couplings which, if allowed, 
would have upset the missing VEV pattern of Eq. (1). Their absence to all orders thus guarantees 
that the pattern of Eq. (1) is absolutely stable (barring of course SUSY breaking at the TeV scale 
which is safe). As far as the color-triplets are concerned, since r^r 7^ in Eq. (6) for the triplets, 
all four pairs become super-heavy, just as desired. 

The two massless Higgs doublets which emerge from Eq. (6) represent the MSSM doublets 
and hd which acquire light masses after SUSY breaking. Let us denote the down type doublets in 
H, H' , C and C by H^, Hj, and Cj respectively, and hkewise the up-type doublets. It is easy 
to see from Eq. (6) that hu is composed entirely of if - i.e. hu — H^, while hd is a mixture of four 
components Hd, Hj , Cd and ■ In particular, the weights of hd in if, ii', C and C are given by 
H D cos 7 • hd, H' D ^y!^ ^ cos 7 ■ hd, C Oi ^gy^' cos 7 • hd and C" D cos 7 • hd- The angle 7 is 
determined in terms of the parameters of the superpotential. It is related to the MSSM parameter 
tan /3 as tan /3 — ^ cos 7. Note that, unhke in many SO{10) models, the MSSM parameter tan /3 is 
not required to be large here. It would turn out that conservative upper limits on proton lifetime 
would correspond to smaller values of tan /3. 

The //-term, the coefficient of h^hd term of MSSM superpotential, is generated within our model 
in a simple way. In the unbroken SUSY limit, //-term is zero since terms such as are forbidden. 
After SUSY breaking, the adjoint ^4(45) develops a VEV~ rususy along its I^r direction, correcting 
the zeros of Eq. (1), which generates the /<-term. This occurs since the inclusion of the soft 
SUSY breaking terms induces VEVs ~ rUsusy foi' the fields C and C' along their z/'^-like scalar 
components. These will trigger the VEV (~ rUsusy) of ^(45) in the I^r direction. From Eq. (4) we 
obtain XiHAH' — > nisusyhuhd, and thus /t ~ mgusy ~ TeV, independent of the integer k. Thus the 
present setup provides a simple and elegant solution to the /t problem without any new ingredients. 

Using Eq. (6), for the four heavy triplets Tj and three heavy doublets Di (coming from four 
pairs of (5 -|- 5)'s of SU (5) in H, H' , C, C , C, C') we derive the following mass relations: 

Md,Md,Md, ^ 9 J_^(m-^\ =^ (7) 

Mt^Mt^Mt^Mt^ 4Mefrcos7 ' M^e v ^ Ai Afa^ ' ^ ' 

y? Q?' 

In our model, for A; = 5 we have Mcfj = ^^^M^. Putting = M\jg^ with ^ 9gvt ~ 
0.63(1 ± 0.10), Mx ~ (0.6 - 1) X 10^^ GeV (see discussion after Eq. (14)), and taking an explicit 
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solution^ for the VEVs, z ~ 0.17M* and s ~ M^^/IQ, together with natural values of the couplings 
Ai (1/4 — \/2) and Xh' ~ (1/4 — 2), we estimate 

Meff ~ (5 X 10^6 - 6 X 10^^) GeV . (8) 

The mass scale M^s will control the d = 5 proton decay amplitude (see e.g. Ref. [9]). It would also 
enter the threshold corrections. Note that Meff, which does not represent the physical mass of any 
particle, can naturally exceed even 10^^ GeV. 

Now, the multiplets A, C, C", C and C' contain three pairs of (10 + 10)'s of SU{5), which get 
masses through Eqs. (2) and (3). Their mass matrix is given by: 



A 



^10 / - \ 
























X2 




Mc 



(9) 



with ^ = {u\q,e^), = (2,1,0), Mvj, = (0,0, Ms/2), where Xi,2 = 4ai,2^c/M^ and Fi,2 are 
defined after Eq. (6). Ms = 2A'^Ma/(6Aa + A'^) yields the mass of the color octet and SU{2)l 
triplet in A(45): Ms = Mg = 2M3. We see from Eq. (9) that two pairs of {u", q, e^)-like states 
are massive. The third massless pair (10 + 10 of SU{5)) is eaten by the corresponding massive 
gauge supcrficlds of SO{10)/ SU{5). Denoting the masses of the components of (10 + 10)'s by curly 
symbols (e.g. = M{u1) etc.), we derive from Eq. (9): U^U^ = 11^2(4+^^), Q1Q2 = 11^2(1+^^), 
S^S^ = YiY2p^ where ^ ^ = ^ = p2 ^ p2 |;L - ^^|. These expressions will be useful 
for the computation of threshold corrections in the model. 

The masses of the heavy gauge boson superfields corresponding to the broken generators of 
50(10) are given by (see e.g., the second paper in Ref. [8] and [9]): M^(X, F) = g'^a'^ = M|-, 
M\X',Y') = Mi(l + p2), M2(K,c,^c) = Mi(4 + p2)^ M^{Vec^e^) = Mip2 where g is the unified 
gauge couphng at the GUT scale. Given the p and p dependence of the masses given above, we 
see that, except for the e'^-like states, threshold corrections from all other states in the (10 + 10) 
matter sector cancel precisely against those from the corresponding states in the gauge sector. 



^Details of estimating the VEVs based on explicit solutions to the Da = condition will be presented in a 
forthcoming longer paper [16]. 
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An accidental = 4 supersymmetry present in the model (the gauge bosons and three pairs of 
matter fields in the (10 + 10) sector form an A" = 4 SUSY gauge multiplet) is responsible for this 
cancelation. This results in an enormous reduction of the parameters, rendering the model very 
predictive for proton decay. 

We have presented the whole spectrum of the theory, except for the singlet sector, which is 
not relevant for the calculation of threshold corrections. We have however, analyzed the singlet 
sector and verified that there are no unwanted pseudo-Goldstone states in the model. While it 
might appear that there is a C/(l) symmetry associated with the "integer" k in Table 1, it turns 
out that this is a linear combination of U{1)a and B — L, and its breaking does not lead to a 
pseudo-Goldstone boson. 

The evolution of the three gauge couplings in the model with momentum is shown in Fig. 1, 
which takes into account all the threshold effects. It is clear from Fig. 1 that the three couphngs 
merge into one at a unification scale Mgut ~ 10^^ GeV. Furthermore, we see that the unified 
5*0 (10) gauge coupling remains perturbative to scales well above Mqut- This is a desirable feature 
which not all SO {10) models have. 

3 Novel correlation between d = b and d = 6 proton decays 

Writing down the three RG equations for Q;]^2,3 and eliminating the unified gauge coupling ao we 
obtain 

(2) 

We have taken GUT scale threshold corrections in one loop approximation. The quantities A • 
include weak scale threshold corrections and 2-loop running effects for the gauge couplings, including 
Yukawa interactions. Their values depend on the SUSY particle spectrum. We carry out our 
analysis within the minimal N — 1 SUGRA scenario [17] with family universal parameters. While 
we vary these parameters to draw our conclusions, for concreteness, we consider the set of values: 
{tan;g,mo,mi/2,/i} = {3,1448.2 GeV, 155.93 GeV, 1 TeV}, which corresponds to rUq ~ 1.5 TeV 
and m^^ ~ 130 GeV . For these values we obtain A- = (0.6093,0.4079, 1.167). Using the third 
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iT I I I I I I I I I I I I I I I 

10^ 10^ lO"* 10^ 10^ 10^ 10^ 10^ 10^°10^M0^^10^h0^''l0^^10^^10^^ 

ix (GeV) 

Figure 1: Evolution of the three standard model gauge couplings in the present SO{10) model 
including threshold corrections. We have used a^i^Mz) = 0.1176 and assumed an mSUGRA spec- 
trum with {tan/3, mo, mi/2, /i} = {3, 1448.2 GeV, 155.93 GeV, 1 TeV} (corresponding to nig = 1.5 
TeV, my^, = 130 GeV), and have taken p = 4, r = 1/250, Meg = 4 ■ 10^^ GeV, 1^1,2 = 2Mx/45 for 
generating this plot. 
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RGE equation for the a^^ we obtain for the present model aciMx) — a^^Mx) — 1/20. The 
parameter k, = Mg/M^ = 2 in our model (as opposed to k = 1 in SUSY SU{b)). 

It is important that the ratio r = M-^/Mx entering into Eq. (11) is constrained by symmetries 
of the model. Using expressions for Ms, Ma and Mx presented earlier, we obtain: 



Mx ~ ^2 M, ^ U5 3007 



(12) 



The range for r is obtained by noting that is allowed by symmetries of the model and thus 
naturally expected to be of order one. We have thus taken, 1/5 < A'^ < 2 (say), and have set 

~ gl^^ ^ 0.63(1 ± 0.10), Mx ~ (0.6 - 1) x 10^^ GeV (see discussion after Eq. (14)), while 
M* ^ Mpi ~ 2 X 10^^ GcV. This restriction on r will be an important ingredient in the derivation 
of an upper limit on r~^(p — )■ e"^7r°). 

Eliminating p/p from Eqs. (10) and (11) we obtain a correlation between Mgs and Mx for a 
given r: 

- (^)i^) {^) m t'-'^i - } • <-) 

where rj^ ~ 0.6 accounts for the running of cos 7, and Sa^^ denotes the deviation of a^^ from 
its central value of 1/0.1176. Note that the curly bracket on the right side of Eq. (13) is fully 
determined for any given choice of the SUSY parameters and a^i^Mz). It turns out to be only mildly 
dependent on variations of mo and mi/2. Since Meg and Mx respectively control d = 5 {p ^ vK^) 
and d = 6 (p — > e"'"7r°) decay amplitudes, Eq. (13) in turn provides a correlation between the 
rates of these two otherwise unrelated decay modes. Such a correlation exists in minimal SUSY 
SU{5) as well [18], but that leads to predictions for a3{Mz) and d — 5 proton decay rate which 
are inconsistent with experiments [19]. Exceptions to this conclusion has been suggested in Ref. 
[20] which uses higher dimensional operators. This however leads to large threshold corrections, 
making the apparent unification of gauge couplings with low energy SUSY somewhat coincidental. 
For review of proton decay in SU{5) and in alternative scenarios see Ref. [21]. 

Now, using expressions for proton decay rates (see below) one finds that the empirical lower 
limit on T^^{p — )■ vK^) requires that M^s > 2.91 ■ 10^*^ GeV (for reasonable scenarios for the 
Yukawa couplings, see discussions in Sec. 4.1), while that on T~^{p e+vr^) requires (owing to 
Eq. (13)) r < 1/150. Using the particular choice of SUSY parameters stated above, and the ranges 
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Figure 2: 



Correlations between Mes 



and Mx for {tan/3, mo, mi/2, /u} 



{3,1448.2 GeV, 155.93 GeV, 1 TeV} (corresponding to = 1.5 TeV, = 130 GeV), 

and a3{Mz) = 0.1176. (a): r = 1/200. (b): r = 1/250. (c): r = 1/300. The vertical and 
horizontal dashed lines correspond to the experimentally allowed lowest values of Mx and Mes 
which arise from limits on r~^(p — t- e+7r°) and T~^{p — )■ vK^) respectively, for central values of 
relevant parameters (see text). 

for Meff and r as given in Eqs. (8) and (12), we therefore illustrate our results on the correlation 
between M^s and Mx in Fig. 2 by confining to the ranges M^s — (2.91 — 6) x 10^^ GeV and 
r ~ (1/200 - 1/300).^ 

Two crucial differences between our model and that of minimal SUSY SU{b) (see Ref. [18]-[21]) 

are: (i) Mh^ of 5'f/(5) is replaced by Meg cos 7 which can be significantly larger than Mqut (for 

tan/3 > 3 say); and (ii) whereas the range for r is severely restricted by symmetry considerations 

for the S'O(IO) model (see Eq. (12)), this is not so for the minimal SUSY SU{^) model. These 

two distinctions make the S'O(IO) model presented here more predictive for proton decay into e^vr^ 

on the one hand, and viable on the other. In particular, as discussed below, with a3(M^) being 

^While such large values of McS lie in their natural ranges, they correspond via Eq. (10) to rather small values of 
p/p ^ 10~^. Thus p/p ^ 10^^ is the only parameter of the model whose smallness remains unexplained on grounds 
of naturalness. 
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consistent with experiments, the d = 5 proton decay rate is in full accord with the experimental 
limit. 



4 Nucleon decay 

There are two main mechanisms for proton decay corresponding to d = 5 and d = 6, which 
respectively yield p DK^ and p c+tt^ as the dominant decay modes. Although apriori these 
two modes are largely independent, owing to the correlation given in Eq. (13) and Fig. 2, they 
get linked in our model such that the observed lower limit on the inverse decay rate of either mode 
implies an upper limit on that of the other. The latter is found to be especially constrained for 
the p — !■ e^Ti^ mode. The rates for (i = 6 decay modes p — )■ e^7r° and p i>7r^, which are largely 
independent of the details of Yukawa couplings and SUSY spectrum, are given by: 

+ ^ + r(p^^vr+)^2r(p^e%°)^^^, (14) 

where f{p) — 4 + (1 + 1/(1 +p^))^. Here an denotes the hadronic matrix element. Recent 
lattice calculation yields au — 0.012 GcV'^ at /i = 2 GeV [22]. D and F are chiral lagrangian 
parameters with D ~ 0.8, F ~ 0.47. gx denotes the effective X, F boson coupling at Mx- 
The correlation curves (see Fig. 2 for a representative case) restrict Mx in the range of about 
(6 — 10) X 10^^ GeV. Taking an average of the three gauge couphngs, which nearly unify at Mx, 
lying in the range as given above (see Fig. 1), we obtain adMx) = gx/i^T^) — (1/20)(1 ± 0.1), 
where the error reflects variations in the GUT scale spectrum or cquivalcntly in the parameters of 
the superpotential lying in a natural range. The function f{p) varies between the limits 8 and 5 
as p varies from to oo; correspondingly one obtains r(p e+7r°)/r(p i/7r+) ~ (1, 1.4, 2.5) for 
(p< 1/3 1 1). The case of p — > oo represents the SU{5) limit. If this branching ratio 

is measured to be significantly smaller than 2.5, that would be strongly suggestive of S'O(IO) (as 
opposed to SU{5)) grand unification. The quantity Ar in Eq. (14) denotes the net renormalization 
of the d = 6 operator, including short {Ag"^ ~ 2.22) and long distance effects {Ri ~ 1.25). In our 
model Ar ~ 2.78. From Eq. (14) we get: 

1/20 V/ Mx y 

aaiMx) ) V6.24 x 10^^ GeV J ' 

(15) 
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Allowing for uncertainty in by ±25% (see discussion in Ref. [22]) and that in adMx) by ±10%, 
and letting p vary in the theoretically favored range of p ~ 1 to p = 10, we see that the empirical 
lower limit on F^^ip e+7r°) > 1.01 x 10^^ yrs [23] requires Mx ^ (6.26 • lO^^GeV) • (1 ± 0.14), 
where for the central value we have used p = 4, and the errors are added in quadratures.^ 

Without further theoretical constraint on Mx, given that r^2g(p — > e"'"7r°) oc M^, there has 
been considerable uncertainty in the literature so far on T~^{p e+vr"), which has been quoted 
to lie in the range ~ 10^^ — 10^^ yrs [24, 20]. Such a range corresponds to (apriori reasonable) 
guesses on Mx ~ Mgut x (1/3 - 3) ~ (0.7 - 6) x 10^^ GeV, the higher values of Mx being allowed 
by letting r be arbitrarily small (^ 10~^ say). Lifetimes much exceeding (few — 10) x 10^^ yrs 
would, however, be inaccessible to next generation proton decay experiments. As discussed below, 
the correlation Eq. (13) together with the restriction on r (see Eq. (12)) would provide a much 
stronger constraint on r~^(p — > e"'"7r°)inax, which is fully accessible. To obtain an upper limit on 
Mx and thus on r~^{p — >■ e"'"7r°) we first need to discuss d — 5 proton decay. 

4.1 Fermion masses and p vK^ decay rate 

In order to investigate d — b proton decay, the Yukawa sector needs to be specified. We have 
found a self-consistent picture for fermion masses and mixings in the present setup in the same 
spirit as in Ref. [25, 9]. Here, for the sake of completeness, we present only the gist of this 
picture, we will return to a more complete presentation in [16]. We introduce a non-Abelian 
flavor symmetry (the quaternionic group) and assign the matter fields 16i,2 as a doublet of 
this group, 1^ = (161,162). I63 transforms trivially under (^4. Two Q4 doublet flavon fields 
^,1^ both with VEVs along (1,0) direction are also utihzed. The Q4 symmetry also enables 
us to successfully address the SUSY FCNC problem [28]. With the U{1)a charge assignments of 
Q(l^) = -(1/2 ± 3/A;), Q(t) — 3/k and Q(^) — 7/k, the relevant operators, in accord with 
the symmetry 5'O(10) x U{1)a x Z2 x Q4, which generate effective Dirac Yukawa couphngs, are: 
163163//, il^lQsH, ^T^T^H, ^Thlhc, ^ {Th ■ 163± I63 • T^) C, and -^(^i^)(l^T^)C". 
The higher order operators, suppressed by powers of l/M^, and in the last two cases by 1/(2')^ as 
well, may be generated by quantum gravity and in part by exchange of additional heavy vector-like 
^This is only to indicate a reasonable range for {Mx)min, which, however, is not used for our explicit predictions. 
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states ( see [16] for details). The resulting mass matrices for the quarks and charged leptons at the 
GUT scale have the form: 

U2 d1{ei) (^2(62) 





( 



U3 



\ 



e' 
-e' 
a 



\ 



m 



di{ei) 
dz{el) 



( 



-Kd(e)e' - rl 



\ 



d\{ez) 


a + Kd(e)e 
1 



(16) 



where m^^ = cos 7/ tan ^5, = 1 and k^. — 3. Eqs. (16) provide a constrained system with 
fewer parameters than observablcs. A consistent fit for all masses and mixing parameters as well 
as observed CP violation is obtained with the choice o = 0.0508, e = —0.0188 + 0.0333i,e = 
0.106 + 0.0754i,e' = 1.56 • 10'^, rj' = -0.00474 + 0.00177i,^^2 = 0.014e^-^* at the GUT scale. 
Following renormalization in going down to low energies (with mt{mt) — 160 GeV and tan^ = 3), 
these values reproduce the central values of the charged lepton masses. In addition, for the quark 
masses we obtain 



m„(2 GeV) = 3.55 MeV, m^{mc) = 1.15 GeV, 

md(2 GeV) = 6.45 MeV, m,(2 GeV) = 137.6MeV, mb{mb) = 4.67 GeV. 
For the CKM mixings we obtain &t n — Mz, 



:i7) 



\Vus\ = 0.225 , \Vcb\ = 0.0414 , iKbl = 0.0034 , \Vtd\ = 0.00878, r} = 0.334 , p = 0.12 . (18) 

Thereby wc get sin 2/3 = 0.663. All these are in a good agreement with experiments. 

Let us now briefly discuss the neutrino sector. The relevant operators, responsible for generating 
heavy Majorana masses for the right-handed neutrinos are: .fLs?! l^^C*^, ^A-^Li ^ ' 1^1636'^, 



and 



I63C . Here, M^, M^i and Mat" represent masses of additional singlets which turn 



out to he in the range of (few - 100)Mgut [16]. We assume that in the first of these couphngs the 
Qi contraction is in the V channel. The heavy Majorana mass matrix Mr is given by: 



R 



6 

6a 
a 1 



Mn 



(19) 
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The Dirac mass matrix at GUT scale can be obtained from M„ (see Eq. (16)) by the replace- 
ment e' — > — 3e'. We can take the two dimensionless parameters (a, 6) and the mass parameter Mq 
as input to fix ^Am^, 612 and ^23- Two observables, viz., yAm^^/Am^ ~ m2/m3 and ^13, will 
then be predictions of the model. The structures given in Eq. (19) are vahd at GUT scale. Apply- 
ing renormalization, including threshold effects due to the different v'^ masses, with 612 — 30° and 
6*23 — 43° as inputs, we obtain 777.2/7713 ~ 0.13 and 6*13 ~ 3.6° as predictions. Such a fit is realized 
by choosing a = 0.0252e-° °i^\ h = 1.61 • lO'^e-^-^^^S and Mq = 1.89- 10^^ GeV. The corresponding 
masses are (M^i, Mr2, Mm) = (3.04 x 10^ 1.2 x 10^°, 1.79 x lO^^) Qgy. These results include all 
the relevant RG running effects. One sees broad, although not precise, agreement with data. We 
consider this fit, which provides large neutrino oscillation angles, together with small quark mixing 
angles as well as observed CP violation as fairly successful and highly nontrivial, especially in a 
quark-lepton unified framework with a stabilized doublet-triplet splitting. 

With the Yukawa couplings specified, the inverse of the sum of partial decay widths, in p — > PX"*", 
is computed to be: 



T^l^ip j?K+) = 3.5 X lO^^yrs 



'0.012GeV^\ /6.91\ /1.25\2/ 



\/3h\ ) \A^s ) \Rl J V3.38 x lO^^GeV 



'cff 



X 



X 



777, \Vl30\ V 3.1 



l.STeV; \m^) \K^^, 



Here i^^=5 denotes a sum of contributions to the total decay rate from the three neutrino fiavors, 
reflecting the dependence of the d = 5 operator on the Yukawa couplings: K^^^ = \ AyJ[^ + + 
1^4,^^ p. Each individual A^,. receives contributions from three types of diagrams leading to the d = h 
operator: (a) those with only the first two families in the external legs, (b) those having the quark 
doublet of the third family in just one external line, and (c) those having the same as in (b) in 
two external lines. The last two contributions incorporate the short distance renormalization of 
the (i = 5 operator that arises through the running of the top quark Yukawa coupling, from the 
GUT scale to the weak scale. Contributions from all three diagrams are found to be important, 
especially for \AvJ^ and \Ai,^\. The net result is that \Ai,J\ ~ 0{1Q~^), \Aj,^\ ~ \Aj,^\ ~ 0{1) and 
K^=^ ~ 3.1 [16]. A'g in Eq. (20) denotes the short distance RGE factor for the d — b operator, 
corresponding to the running from the GUT scale to the weak scale, that arises purely from the 
gauge interactions, without the effects of the top quark Yukawa coupling. Note that A'g defined 
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here differs from tlie RGE factor As defined conventionally [18] in that As includes the effect of 
the running of mciris in going from low energies to the GUT scale, Ag does not. Thus, As — AgJ, 
where J = {'mc'ms)G\]T / {'fnc'nis) ^ ~ O(10~^) for // = 2 GeV (with low tan^ ~ 3 to 10). 

We can now discuss the derivation of an upper limit on Mx and thus on T~^{p — > e^i:^). 
Owing to Eq. (13) this would correspond to the minimum allowed value of Mcfj. Now, taking 
conservatively rUq < 1.5 TeV and the experimental lower limit > 125 GeV, the observed lower 
limit on V-'^{p VK+) > 2.8 x 10^^ yrs [23] yields (via Eq. (20)): (Meff)min ^ 2.91 x 10^^ GeV. 
This in turn yields by using Eq. (13) (with \/3h\ = 0.012 GeV^ the lowest value of r = 1/300 and 
tan/3 = 3): (Mx)max ^ (5.16, 7.02, 9.45) x 10^^ GeV for a^iMz) = (0.1156, 0.1176, 0.1196). Thus, if 
wc use central values of 03 (M^), an, Ph and ^^(Mx) with p = 4, the upper limit on T~^{p — )■ e+7r°) 
(using Eq. (14)) would be 1.61 x 10^"^ yrs. If the uncertainties in all these parameters are stretched to 
their extremes, each in a direction so as to prolong the lifetime, the stated upper limit could increase 
by atmost a factor of 10.8. Considering that all the uncertainties having such extreme values, and 
in the same direction, to be very unlikely, we would regard something like the geometric mean of 
the two upper limits corresponding to the central and extreme values of the parameters to be a 
more realistic, yet conservative, upper limit for the lifetime. We thus predict: 

T^Up ^ e+7r°) < 5.3 x 10^^ yrs . (21) 

If ruq < 1.5 TeV, or m^^ > 125 GeV, or r > 1/300, or tan^ > 3, the upper limit would of course 

decrease further significantly.^ Thus, the upper hmit shown above on r~^(p — > e+Tr'^), stemming 

from Eq. (13), is a robust and novel feature of the model. The predicted lifetime is accessible to 

proposed megaton size water Cherenkov (or equivalent) detectors. 

Reversing the procedure given above, we can derive an upper limit on Mes and thereby on 

T~^{p — > VK'^). Owing to Eq. (13), this would correspond to the minimum allowed value of Mx and 

r. Using central values of {a^l and adMx) with p ^ 4 (for concreteness) , the observed lower limit 

onr-^(p ^ e+7r°) > l.OlxlO^Vrs [23] yields via Eq. (15): (Mx)min > 6.26x10^^ GeV. This in turn 

*While r~^{p — >• e+TT^) given by Eq. (14) does not explicitly depend on m,, m^, r and tan/3, the upper limit 
on Mx and thereby T~^{p — >• c+tt") does depend on these parameters via the correlation Eq. (13). The latter 
relates (Mx)max to (Meff)min and thereby to the empirical lower limit on T~^{p — >■ 17K~^) which depends on rriq 
and niyy. Because of this, the upper limit given in Eq. (21) should in fact be multiplied by an approximate factor 
(mg/1.5 TeV)8/3 (125 GeV/m^)^/^ [(l/300)/r]4/3 (3/tan^)4/3. 
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yields by using Eq. (13) (with = 1.5 TcV, = 130 GeV and the lowest values of r ~ 1/300): 
(Meff)max < (1-627, 4.105, 10.02) X 10^9 GeV(3/tan/3) for asiMz) = (0.1156, 0.1176, 0.1196).^ 
If we use central values of the parameters and the spectrum as noted above, with p ^ A and 
tanP > 3, the upper hmit on r~^(p VK'^) (using Eq. (20)) would be 5.16 x 10^^ yrs. Allowing 
for uncertainties in the parameters in a combined manner, analogous to the case of o? = 6 lifetime, 
we thus obtain^" 

V-\p ^ VK^) < (3.1 X 10^^ yrs) x {j;^)' {^^^^^ (S/tan/?)^ . (22) 

In Eq. (22) the mild dependence of the curly bracket of Eq. (13) on rriq and m^i^ is not exhibited. 
The actual lifetime is hkely to be significantly lower than few xlO^^ yrs if Mes is not stretched 
to its upper limit (corresponding to e.g., Mx > (Mx)mm, or r > 1/300, or tan/5 > 3 and/or 
0^3(^2) < 0.1196), or if rrig < 1.5 TeV, or m,y > 130 GeV. We thus find that not only the 
p — >■ e+7r° mode, but very likely even the p — >■ i/K^ mode should be observable by improving the 
current experimental sensitivity by about a factor of five to ten. 

Some important details concerning the present work, including those pertaining to the issues 
of fermion masses and mixings, and some variants as regards the cancelation of the U{1)a Fayct- 
Iliopoulos term, will be presented in a forthcoming longer paper [16]. 

In summary, we have presented a class of supersymmetric 5*0(10) models with low dimensional 
Higgs system that fully resolves all the naturalness issues of doublet-triplet sphtting, including 
stability against higher order operators, generation of //-term of order nisusy, and proton stability. 
The threshold corrections in these models arc found to depend only on a few effective parame- 
ters, making the scenario very predictive. An intriguing feature of these models is the correlation 
equation and the corresponding constrained upper limit on r~^(p e"'"7r°). We find that in this 
class of models proton decays into both e'^n^ and very hkely PK'^ as well should show with an 
improvement in the current sensitivity by about a factor of five to ten. The building of a megaton 
water Cherenkov detector (or equivalent) would thus be most welcome. 

We would like to thank Takaaki Kajita and Edward Kearns for helpful communications. K.S.B. 
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suppressed, owing to symmetries present in the model, so as not to disturb the upper limit quoted in Eq. (22). 
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